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Precision laser spectroscopjS' of cold and trapped molecular ions is a powerful tool for fun¬ 
damental physics, including the determination of fundamental constant^, the laboratory 
test for their possible variatioi®, and the search for a possible electric dipole moment of 
the electron^. While the complexity of molecular structure facilitates these applications, 
the absence of cycling transitions poses a challenge for direct laser cooling, quantum state 
controP®!, and detection. Previously employed state detection techniques based on photo- 
dissociatioiJi^ or chemical reactioni^ are destructive and therefore inefficient, restricting 
the achievable resolution in laser spectroscopy. Here we experimentally demonstrate non¬ 
destructive state detection of a single trapped molecular ion through its strong Coulomb cou¬ 
pling to a well-controlled co-trapped atomic ion. An algorithm based on a state-dependent 
optical dipole forc^ (ODF) changes the internal state of the atom conditioned on the inter¬ 
nal state of the molecule. We show that individual quantum states in the molecular ion can 
be distinguished by their coupling strength to the ODF and observe black-body radiation- 
induced quantum jumps between rotational states of a single molecular ion. Using the de¬ 
tuning dependence of the state detection signal, we implement a variant of quantum logic 
spectroscopj23l^ of a molecular resonance. The state detection technique we demonstrate is 
applicable to a wide range of molecular ions, enabling further applications in state-controlled 
quantum chemistr^ and spectroscopic investigations of molecules serving as probes for in¬ 
terstellar cloudi^i^l^ 

One of the salient features of trapped ion systems is that the universal Coulomb interaetion al¬ 
lows strong eoupling of diverse quantum objeets, sueh as different speeies of atomie ions or atomie 
and moleeular ions. Being able to perform quantum logie operations e.g. in the form of gate^^H^Hl 
between the quantum objeets has proven a powerful tool for quantum information proeessing and 
quantum simulations in sueh systems. It also allows eombining the advantages of different atomie 
speeies. Quantum logie speetroseopy is one sueh applieation in whieh the high degree of eontrol 
aehieved over seleeted atomie ions is extended to speeies over whieh sueh eontrol is laokin^f^HI. 
Here, we demonstrate for the first time quantum logie operations between a single moleeular ion 
and a eo-trapped atomie ion, making a wide range of moleeular ions aeeessible to this highly- 
developed toolbox. The presented teehnique allows the investigation of single moleeules in a well 
isolated system avoiding disturbanee from the environment, whieh is the limiting faetor in other 
implementations of single moleeule speetroseopy sueh as surfaee enhaneed Raman speetroseopy 
(SERSjSH. 
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Quantum logic operations between atoms are based on state dependent forees often indueed by 
laser fields. The same approaeh is applieable to moleeular ions. The eoupling is now distributed 
over many ro-vibrational transitions, whieh weakens the interaetion strength for individual transi¬ 
tions (see Methods). The energy splitting between quantum states allows the ODF to distinguish 
between them, faeilating rotational state deteetion. By tuning the frequeney of a laser field near 
one of the moleeular resonanees, an ODF is exerted on the moleeular ion if it is in one of the states 
eoupled by the eorresponding transition. If the strength of this foree is oseillating with a frequeney 
elose to one of the seeular motional frequeneies of the two-ion erystal in the trap, it is resonantly 
enhaneed, generating eoherent states of motion, whieh ean be deteeted on a eo-trapped atomie 
iorP^. 



Figure 1 : Coupling strength of an optical dipole force to the atomic and molecular ions, (a) 

Partial level seheme of the relevant atomie and moleeular states and their laser eoupling to the 
lattiee laser. The inset shows the moleeular substrueture for the states involved in the sequenee. 
The number on the transition is the experimentally implemented geometrie eoupling eoeffieient. 
(b) The two lattiee beams form a moving interferenee pattern leading to a temporally and spatially 
varying light foree on both ions, (c) Rotational state (J, mj = 0) and detuning-dependent eoupling 
strength of the dipole foree interaeting with a moleeular ion on the P-braneh ( J —)■ J — 1, 

solid line) and R-braneh (J —)■ J -I- 1, dashed line) of the —)■ eleetronie transition. 

The blue dotted line is the eoupling strength to the ion. 


We extend the simple deteetion seheme deseribed above to deteet the small ODF arising from 
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Figure 2: Molecular state mapping to the motional qubit, (a) Circuit description of the mapping 
process. For ideally chosen parameters AMgH and 50, as shown, the circuit is equivalent to a CNOT. 
The action of the ODF on the motional qubit is represented by the torque Fqdf = 

(b) For J ^ 1 only the atom interacts with the lattice. The ODF creates a torque Fug on the Bloch 
sphere (purple arrow), (c) If J = 1, the additional force on the molecule (torque FMgH depicted in 
blue) is added to the force on the atom leading to a modified population transfer. 
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a laser field tuned near the X^E+(J =!)—)■ A^E+(J = 0) resonanee (detuning AMgn) in a 
moleeular ion in the presenee of a nearly detuning-independent baekground foree from a 
eo-trapped atomie ion (Fig. [^). Blaekbody radiation (BBR) induees transitions between 

rotational states J, resulting in a thermal rotational state population with a maximum at J = 4 for 
temperatures near 300 K. We have ehosen to probe the J = 1 state, with a population of around 
8 %, for its long average dwell time of around 3 s (see Extended Data Fig. [^. Vibrational popula¬ 
tion with vibrational quantum number z/ > 0 in is negligible at room temperature owing 

to the large vibrational spaein^^. Applying the ODE resonant with either the IP or OP mode imple¬ 
ments a displaeement operator resulting in eoherent states of motioiJ^, whieh are impraetieal for 
state diserimination due to their non-orthogonality. Instead, we map the moleeule’s internal state 
information to an engineered motional 2-level quantum system (motional qubit) that is eoupled 
by the ODE. This approaeh allows us to suppress the baekground foree from and enables 

advaneed eoherent population transfer sehemes, sueh as eomposite pulsed. 

The state deteetion algorithm starts with both axial motional modes [in-phase (IP) and out- 
of-phase (OP)] of the two-ion erystal eonsisting of a and a ^^MgH^ ion initialized in the 

motional ground state |0)j^ = |0)jp |0 )qp via the atomie iorP^. Two long-lived hyperfine states 
in serve as the atomie qubit denoted by ID^g ^ It)Mg initialized to |t)]y[g. The 

motional qubit is prepared in 1),)^^ = |l)ip |0 )qp by driving a 7r-pulse on the blue sideband tran¬ 
sition (BSB) of the IP mode on the atomie ion while ehanging its internal state from |t)Mg to 
It)Mg- The ODE is implemented by interfering two eounterpropagating laser beams along the ax¬ 
ial direetion of the ion trap, forming an oseillating ID optieal lattiee. By adjusting the relative 
detuning 6 between the two beams to the differenee between axial IP and OP mode motional fre- 
queneies, 6 = wop — cuip = Ci;ip(\/3 — 1), a Raman eoupling involving only the motional states 
lt)m lt)m = |0)ip |l)op implemented, resulting in Rabi flopping between the two eomposite 
states (Fig. [^and Extended Data Fig. [^. The interaetion between the lattiee laser field and the 
motional qubit is deseribed by the Hamiltonian (see Methods) 

H = hQesViPVop (^flipflop + OipOopj ) (1) 

with the ereation operator a|p(op) IP(OP) mode. The effeetive Rabi frequeney is given by 

^eff = ^Mg-y/l — 2i?(AMgH) COs{6(j)) + i?(AMgH)^, (2) 

where R{AMgH) = OMgH(AMgH)/^Mg is the ratio between the Rabi frequeneies of the Raman 
eouplings in the moleeular and atomie ion, depending on the detuning AMgH- The distanee between 
the ions determines the relative phase 50 between the forees on the ions, and r/ip^op are the Eamb- 
Dieke parameters for the OP and IP modes, respeetivel)^^. In a Bloeh sphere pieture, the effeet 
of the ODE on the atom and the moleeule ean be illustrated as the veetor sum of two torques 
-^ODF = -^Mg + -^MgH aeting on the motional qubit (Fig.[^,C). In the experiment, the time fiatt in 
whieh the lattice laser interacts with the ions is adjusted such that for DMgH — 0 (i.e. the molecular 
ion is not in J = 1) a tt- pulse is performed on the motional qubit (Fig.[^). If the molecule is in the 
probed state (DMgH 7 ^ 0), the additional force changes the population transfer (Fig. [^). This step 
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transfers the internal information of the molecular state to the motional qubit. The second BSB 
TT-pulse maps the motional state population (containing the molecule’s internal state information) 
to the atomic qubit. The population in state follows 


■^IDivig 9 


1 + COS f vr^l - 2i?(AMgH) cos(50) + /((Angn)' 


(3) 


and is detected using state-dependent fluorescence on the ion, where fluoresces and 

It)Mg The full sequence is shown in Extended Data Fig.j^ 

The calculated signal is shown in Fig. [3h for different phases 6(j) and detunings AngH- 
The maxima correspond to rotations of 27 m on the Bloch sphere of the motional qubit with n E 
IN. In particular, at the outermost peaks (see Fig. |^) the forces from the atomic and molecular 
ion cancel or add up to 27r exactly, depending on the sign of AMgH- There is a trade-off in the 
choice of lattice laser detuning from the molecular resonance: if the detuning is too large, the 
ODF will not be detectably different for a molecule in the rotational state of interest compared to 
a molecule in any other state. For a detuning too small, spontaneous Raman scattering will change 
the rotational state before it can be detected. Under ideal experimental conditions, it is possible 
to choose a combination of the phase 50 and the detuning AngH to perform state detection in a 
single experimental repetition (cycle), since it is free of quantum projection noise = 1 with 

the molecules in J = 1 and = 0 if it is not). In this case the experimental sequence can 

be understood in terms of a CNOT gate acting on the motional qubit with the molecular rotational 
qudit I J)MgH as a control (Fig.|^): 


^ l)MgH 'l')m 

7^ 3)MgH t)m ’ 

(4) 

^ l)MgH t)m 

0)jn ; 

(5) 

~ 3)MgH i)m 

~ 3)MgH i)m ; 

(6) 

~ 3)MgH t)m 

~ 3)MgH t)m • 

(7) 


A typical measurement with positive detection event can be seen in the inset of Fig. [^. If the 
molecular ion is not in the probed rotational state, the vr-pulse on the motional qubit is successful 
and we observe low fluorescence (Pij.)^^ ~ 0). After entering the J = 1 state the additional force 
on the molecule will change the vr-time on the motional qubit and we observe fluorescence on 
the ion (Pq)^^ > 0), until the rotational state changes again. The signal corresponds to a 

BBR-induced quantum jump of the rotational state in the molecular ion. The signal background 
of Pq)^g ~ 0.2 and reduced contrast arises from experimental imperfections in implementing the 
logic gates (see Methods). Since the algorithm implements a quantum non-demolition measure¬ 
ment of the molecular state, we reduce quantum projection noise by averaging each data point over 
30 experimental cycles with a duration of ~ 9.3 ms per cycle. An experimental event was identi¬ 
fied as a positive detection event if at least three consecutive data points were significantly (1.5 a) 
above the background (see Methods). In a total measuring time of approximately 23 h we have 
observed 18 positive detection signals of up to several seconds duration. The estimated average 
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Figure 3: Non-destructive state detection, (a) Theoretically calculated signal on the 

atomic ion for in the | J = 1, mj = state, and its dependence on the relative phase 

5(f) of the optical force between the atomic and molecular ion and the detuning from molecular 
resonance, (b) Cross section of a, for a phase 5(p 0.6 x vr, corresponding to a trap frequency of 

uyp ~ 2.17 MHz. (C) and (d) show the height of the observed signals as a function of the detun¬ 
ing for the regions indicated in b. A typical measurement result for a black-body radiation driven 
quantum jump into and out of the J = 1 rotational state is shown in the inset. The error bars (68 % 
confidence interval) are dominated by quantum projection noise. The thickness of the theory bands 
accounts for small changes in trap frequency (and thus 5(f)) for the different data points. The grey 
data points are excluded from the analysis in Fig. since they could not be assigned to a specific 
slope of the theoretical prediction, thus preventing unambiguous inversion of Eq. Q. 
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number of scattered photons from the optical lattice beams interacting with the molecular ion dur¬ 
ing a single detection cycle at AMgH = 25 GHz is around 0.04 (see Methods). If we assume that 
a single scattered photon removes the molecule from its rotational state, this results in an average 
single-cycle detection efficiency of 96 %, neglecting imperfections in any of the qubit and gate op¬ 
erations, thus representing an upper bound. Estimates of the expected number of events using this 
detection efficiency together with the mean population and the required averaging time of ~ 1 s 
for a positive signal event agree well with the experimental results (see Methods). 

The lattice light introduces an AC Stark shift on the order of several hundred kHz, which is 
much larger than all magnetic couplings in the molecular ion. Therefore, the substates are labeled 
with mj in a basis with the quantization axis along the electric field direction of the lattice light. 
In this basis, only vr-polarized light is present, coupling only to states with mj = 0 (see inset 
Fig. [T^). The resulting signal as a function of the detuning and molecular quantum state is shown 
in Fig.|^ and d for Scj) ~ 0.6 x vr, corresponding to a normal mode frequency of uyp ~ 2.17 MHz. 
Experimental data for 50 0 is shown in Extended Data Fig.[^ 

The dependence of /((Angn) on the detuning AngH from a molecular resonance allows us to 
perform spectroscopy, since forces arising from coupling to other rotational states are negligible 
over a large range of detunings (Fig. [^). For this, we have combined data taken for different 
detunings AngH and phases 50 by inverting equation Q numerically and accounting for the re¬ 
duced contrast (see Extended Data Fig.[^. The measured Rabi frequency ratios /((Angn) for the 
X^S'''(J = 1) —A^S’''(J = 0) transition are summarized in Fig. together with a fit to the 
wings of a simplified model Rhtif) = A/\f — fo\ with the free parameters A and /o, where we as¬ 
sumed a detuning-independent background coupling to the far detuned transition in the ion 

(Fig.[^). We determine the transition frequency to be /o = 1067.74789(40) THz (68 % confidence 
interval) which is in agreement with a previous measurement of /o = 1067.7473(15) TH^. The 
fitted characteristic detuning for the Rabi frequency ratio is A = 45.98(73) GHz (68 % confidence 
interval). 

We estimate that non-destructive detection of the internal state of a molecular ion with near 
100 % efficiency could take less than 10 ms with a larger detuning of AjvigH ~ 100 GHz and im¬ 
provements in the atomic ion’s state preparation and manipulation (see Methods). This represents 
a significant advance over previous destructive detection techniques requiring frequent reloading. 
The technique will unfold its full potential in high precision spectroscopy of narrow transitions 
using an independent spectroscopy laser. However, while in the present work black-body radi¬ 
ation probabilistically populates the detected state, precision spectroscopy will require efficient 
state preparation schemed^^ll^ or ro-vibrational cooling techniqued^®^. A combination of these 
powerful tools will enable the realization of optical clocks based on molecular ions approaching 
the 10“^® leveP^, where the underlying clock transitions or a combination of transitions can be 
sensitive to variations of fundamental constant^, an electron electric dipole moment (eEDM)P or 
parity violation in chiral molecules. 
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Figure 4: Quantum logic spectroscopy. The Rabi frequeney ratio R = f2MgH(AMgH)/^Mg as 
a funetion of detuning AMgH from a molecular resonance for all data for which Eq. could be 
inverted unambiguously (see Fig. |^,d). Measurement data recorded at ccip 2.17MHz (<50 
0.6 X TT, 9 Points, see Fig. and ccip ~ 2.21 MHz (6(p ~ 0, 7 Points, see Extended Data Fig. 
are indicated by blue circles and red diamonds, respectively. The solid line is a fit to the combined 
data set resulting in /q = 1067.74789(40) THz. The error bars correspond to the errors shown in 
Fig. and Extended Data Fig. [^(68 % confidence intervals). 
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Methods 

Experimental setup The experiments were performed in a linear Paul trap, with typical ax¬ 
ial trap frequency of 2.21 MHz for ^^Mg'*'. We use a hyperfine qubit consisting of the states 
l'l')Mg = 1-^ = 3, mi? = 3) and |t)Mg = 1-^ = 2,mi? = 2) of the ^Si /2 electronic ground state. 
The Raman laser system for coherent manipulation of the ^^Mg"*" qubit via the ^P 3/2 Eevel is 
based on a frequency-quadrupled fiber lasei^. The detuning from the atomic resonance is fixed 
at 9.2 GHz, limiting ground state cooling efficienc)^^ and qubit operations due to spontaneous 
Raman scatterin^^. The same laser is also used for internal state discrimination by applying an 
optical sideband tuned near resonance of the ID^g PP 3 / 2 , F = A, mp = 4) component of the 
D2 transition. A similar laser system addressing the ^Si /2 —)■ ^Pi /2 (Dl) transition is used for 
optical repumping and quenching during sideband coolin^^. 

The light fields creating the optical lattice for molecular state detection and spectroscopy are 
provided by a frequency-doubled dye laser whose frequency is tuned close to the X^E+(J = 
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1) —)■ A^S+( J = 0) transition of with a variable detuning | Angnl ^ 100 GHz from the 

moleeular transition and detuned by AngH ~ 1-5 THz from the atomie D1 transition of 25Mg+. 
Two eounterpropagating laser beams derived from this laser are linearly polarized in the horizontal 
plane with a relative detuning S. They are aligned along the trap axis forming an angle of tt/4 with 
the atomie quantization axis given by a statie applied magnetie field of ~ 0.58 mT. The detuning 
between the two laser fields results in a moving optieal lattiee, providing a temporally and spatially 
varying foree on the two ions with a relative phase 5(f), that depends on the distanee between the two 
ions. Two optieal fibers suitable for UV single mode transmissioEpl are used to reduee the effeet 
of beam pointing fluetuations on the eoupling strength and aid in the initial alignment of the beams 
onto the moleeular ion. The frequeney of the speetroseopy light is monitored by a eommereial 
wavemeter (High Finesse WS/7). We assign an uneertainty of less than 100 MHz to the measured 
value. 


To reduee quantum projeetion noise due to imperfeet single qubit operations, we average over 
30 measurement eyeles to get a single measurement point. This proeedure takes 280 ms. If at 
least three eonseeutive points are signifieantly above the threshold (1.5cr), we assign it to a positive 
deteetion signal from the moleeule. For the determination of the signal height, the first and last 
measurement point were only taken into aeeount if they were at least 2a above the baekground. 
This proeedure improves the signal quality, sinee data points are removed whieh eontain only a 
fraetion of positive deteetion events during the 30 repetitions. 


Hamiltonian for lattice coupling The Hamiltonian for N ions eoupled by a ID moving optieal 
lattiee in the interaetion pieture ean be written as 


H 


N 

E 


exp 


■ N 

E 


IVjk 


a, 


+ h.c. 


( 8 ) 


where j labels the ions and k labels the axial motional modes with frequeney Uk- rjjk is the Lamb- 
Dieke parameter and a|, is the ereation operator, flj = flijfl2j/{2A) is the two-photon Rabi 
frequeney for the two lattiee beams eoupling to an eleotronieally exeited state with detuning A. 
Considering the ease of two ions j G {1, 2} we have to take two axial motional modes into aeeount 
(k G {IP, OP}) denoting the in-phase and out-of-phase motion of the two ions, respeetively. In the 
following, the relative detuning of the two lattiee beams, 5 is ehosen to be 5 = cuop — i^ip- After the 
rotating wave approximation only the seeond order mixed term in the Lamb-Dieke approximation 
survives, resulting in the following interaetion Hamiltonian, mediating the eoupling between the 
two modes: 

H = ~ ^ (fGje''?^^''7jjp?7j^opaipaop + h.c. (9) 

i=i,2 

In adding up the forees on the two ions we have to eonsider the different effeet on both modes and 
take the signs in the matrix for transformation between normal modes and laboratory frame into 
aeeounP^. We find for two ions with equal mass the Lamb-Dieke faetors 


— "72,IP = ^IP5 

= —'72,OP = VOP- 


7i,iP 

7i,op 
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( 10 ) 

( 11 ) 




and the Hamiltonian reads 

H = —hT]ipr]opaipaQp — ^ 26 “^^] + h.c.. (12) 

It is convenient to define 


6(j) = ((> 1 - 02 , (13) 

<I> = 01 + 02. (14) 

The Hamiltonian can then be written as 

H = —e^^^^hrjiprjopaipaQp — 02^”^'^'^^^] + h.c.. (15) 

We introduce the modified Rabi frequency 

Q = , (16) 

with amplitude 

Heff = ^2^1 — 2i?cos(50) + (17) 


where i? = H 1 /H 2 is the ratio between the Rabi frequencies and <l>efr is an effective phase. Identi¬ 
fying f2i = flivigH and H 2 = ^Mg with the Raman Rabi frequencies for the molecular and atomic 
ion, respectively, we arrive at Eq. (2) of the main text. For very large detunings from either the 
atomic or the molecular ion, the effective Rabi frequency reduces to flngH and flivig, respectively. 
The final interaction Hamiltonian is then given by 

H = hQesViPVop (^e'^"®aipaop + e"‘'*’"®aJpaop j , (18) 

which is a Hamiltonian that describes a system consisting of n + 1 coupled levels, where n is the 
number of initial excitations. In the case n = 1 we get a two-level system consisting of the states 
= |l)jp |0 )qp, = |0)jp |l)op which are coupled with the Rabi frequency iles- The action 

of this Hamiltonian can be described as Rabi flopping between motional states. Results from an 
experimental implementation are shown in Extended Data Fig.|^for a two-ion crystal consisting of 
a ^®Mg+ ion and a molecular ion in a rotational state not addressed by the lattice.The reduced initial 
contrast and its decay are attributed to imperfect ground state cooling, off-resonant scattering of 
the Raman beams interacting with the atom, and dephasing of the lattice beams. Starting from both 
axial modes near the ground state, a quantum of motion is added using a BSB 7 r-pulse. The lattice 
coupling is then switched on for a variable time and the motional information is mapped back onto 
the electronic qubit by a second BSB 7 r-pulse (see Extended Data Fig.|^. In such a sequence, we 
are not sensitive to the phase factors which can therefore be neglected in the Hamiltonian, 

resulting in Eq.Q of the main text. Note that an additional phase of n is introduced if the sign of 
the detuning differs for the atomic and molecular resonances, so the relative phase is given by 


(50 


27rdAk 
27i{dAk -f 3) 


(same sign) 
(opposite sign). 


(19) 


where d is the distance between ions and Ak = fci — ^2 ~ 2ki is the difference between the 
wavenumbers of the two lattice beams. 
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Extended Data Figure 5: Rabi flopping between motional states. Implementation of the sequenee 
shown in Fig. for flngH = 0. The duration of the applied optieal lattiee is varied to induee Rabi 
flopping between the motional qubit states and The error bars show the 95%-eonfidenee 
interval of the photon distribution fiP^. 


Coupling of the optical lattice to the molecular ion The coupling strength of the lattice laser to 
the molecular ion can be described by the transition matrix element between two states A"u” J”m'j) 
and \^'K'v' J'm'j) with A as the projection of the angular momentum on the molecular axis, which 
is A' = A" = 0 for both involved states, since we have at a S —)■ S transition without electron 
angular momentum. The electronic part of the wavefunctions is given by = A and = X 
in spectroscopic notation. The vibrational and rotational quantum numbers are given by v and J, 
respectively, and mj is the projection of J onto the quantization axis. Note that the quantization 
axis for the molecule is given by the electric field of the lattice light, since the AC Stark shift causes 
the dominant level splitting. The transition dipole matrix element is usually approximated by 


(ev'jw^A'i dg ie'v"j"m"A") = 




( 20 ) 


where /xxa is the electronic transition dipole moment. The Franck-Condon factor describing the 
overlap of the wave functions of the two vibrational ground states can be calculated from the data 
given in Aymar et al.^ 




FC/ 


/ = 0, v" = 0) = 


n 2 




0.0919. 


( 21 ) 


The Honl-Fondon factor describing the dependence of the transition strength on the rotational 
quantum numbers is given by Hansson and WatsoiP^l 

A“^(J', A' = 0; J", A" = 0) = max(J', J"). (22) 
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The geometric factor describing the addition of total angular momentum of the molecule and the 
angular momentum of the photon is given by 




J' 1 J"\ 
-m'j q 


(23) 


with the Wigner 3j-symbol and the polarization components e'^. Since the quantization axis is given 
by the lattice, it is convenient to calculate the coupling in a basis parallel to the lattice electric field, 
where the polarization components are e~^ = = 0 and e° = 1. 


The Rabi frequency from the interaction of the optical lattice with the molecular ion (Eq. [T8| ) 
starting from a specific rotational state v” J”m'j) is given bjJ23 




MgH 


E 1 E 2 

2h‘^ 2^2^ 

m'j q 
E1E2 fixA qFC 




A 


MgH 


2h^ A 


MgH 


^FC^HL 

EE ^geom 


E\E2 HxA cFC cHL cgeom 
^ 

2h? A MgH 


(24) 

(25) 

(26) 


where Ei and E 2 are the electric field amplitudes of the two counterpropagating laser beams 
forming the optical lattice and AMgH is the detuning from the molecular resonance. Using a 
similar approach for the atomic ion, we get a Rabi frequency ratio that is independent of the 
applied electric field. From the well-known transition strengths in and the fitted value for 

A = 45.98(73) GHz (68 % confidence interval), we extract a dipole matrix element for the molecu¬ 
lar transition of i/xa — 1.2510(99) eoo, which is smaller than the matrix element = 1-779 eao 
extracted from the theoretical data presented in Aymar et We use the experimental dipole 
matrix element for the curves in Fig. and Extended Data Fig. 


Scattering rate and estimated single-cycle detection efficiency We assume that scattering of a 
single photon on the optical transition of the molecular ion will change its ro-vibrational state, 
due to the non-vanishing Franck-Condon factors between different vibrational states. Therefore, a 
small scattering rate is desired. In the following we use a simplified model in which the molecular 
ion is treated as a two-level system to estimate the order of magnitude for the spontaneous emission 
rate of the molecular ion subject to the ODE The number of scattered photons from a two-level 
system with excited state linewidth T 0 interacting for a time It, with an optical lattice formed by 
two laser beams of equal Rabi frequency U and detuned by AMg from resonance is given by 


V f = 

^ SC^TT 


ToU^ 

2A2 

^^MgH 


(27) 


For the implemented quantum algorithm the time fjr is chosen such that the coupling to the atomic 
ion drives a vr pulse between the two motional states, = vr/ {^MgViPVop) ■ Using R = UMgu/^^Mg 
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and l^MgH = (2AMgH), we obtain 


T f = 

■L Srt'TT 


where T o for the molecular ion is given by 


rnvri? 


AngH^ip^oP ’ 


rn = 




3Tieoh, 


/^XA- 


(28) 


(29) 


with the scaled transition dipole moment /ixA given above. For Poissonian photon distribution the 
probability to scatter no photon follows 


Pscatt(O) = e 


(30) 


For a detuning of ~ 25 GHz we expect to scatter on average 0.04 photons, which corresponds 
to a fundamentally limited single cycle detection efficiency of E^s = F’scatt(O) ~ 96 %. For an 
increased detuning of —100 GHz the number of scattered photons on average would be reduced to 
0.0026 for a vr-pulse with the optical lattice. However, to get a sufficiently large signal contrast for 
single-cycle detection it would be necessary to perform a 27r-pulse with the lattice, which leads to 
a single-cycle detection efficiency of 99.5 %. 


The number of positive detection events iVev for a total number iVtotai of trials can be estimated 
from 

Aev = (F;ss)^^^^^‘^Ptherm^t<3A'total, (31) 

where Ptherm is the probability for the molecule to be in the probed state for a thermal distribution 
at 300 K and Pt >3 is the probability for the molecular ion to stay at least for the time of three 
measurement points in the probed state, which was a chosen prerequisite for positive detection in 
the experiment. In our case, we have Ntotai ~ 2 x 10®, Ptherm = 0.08, Pt <3 = 0.8 (see Monte Carlo 
simulation in Extended Data Fig.|^ and Ncyde = 30. For a single shot detection efficiency on the 
order of Pss ~ 95 %, this results in an expected number of points belonging to a detection signal 
of iVev ~ 130, which is in agreement with the observed number of points iVg^P 100 assigned to 
the 18 positive detection events. 

Uncertainties of the fitted parameters We extracted the center frequency and the line strength 
from a weighted non-linear least square fit. The assigned uncertainties correspond to the 68 % 
confidence bounds C, that have been calculated assuming a normal distribution of errors and ac¬ 
cording to C = b± ty/S, where b is the fitted coefficient, t is the inverse of Student’s cumulative 
distribution function for the required confidence and S = (X^X)“^s^, with the mean squared 
error and the Jacobian X of the fitted values with respect to the coefficients. 


16 






Rotational state J 


Extended Data Figure 6: Single trajectory from a Monte Carlo simulation of molecular dy¬ 
namics. (a) The ion initially prepared in the rotational ground state is transfered to higher ro¬ 
tational states due to eoupling to BBR at 300 K. (b) The probability to find the ion in a eertain 
rotational state in the simulation (red bars) follows a thermal distribution. The blue bars are ealeu- 
lated values from a master equation approaeh. The deviation is due to the finite time interval of the 
Monte Carlo wave funetion simulation. (C) The dwell time deereases for higher rotational states. 
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Extended Data Figure 7: Full experimental sequence, (a) Circuit description of the sequence, 
(i) A BSB vr-pulse initializes the atom in the state |t)Mg motional state in the qubit state 

(ii) The ODE rotates the motional qubit controlled by the internal state of the molecule (see 
Fig. 1^. (iii) A second BSB vr-pulse maps the motional state (that contains the information about 
the internal state of the molecule) to the atomic qubit, (iv) The atomic qubit is read out by state 
dependent fluorescence, (b) Pictorial representation of the laser couplings in a simplified level 
scheme. 
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Detuning AmsH (GHz) 


Extended Data Figure 8: Raw data for <50 « 0. (a) Theoretieally predieted signal for 50 0, 

eorresponding to cuip ~ 2.21 MHz. (b) Zoom into the region shown in a with the measured 7 data 
points included as red diamonds in Fig.|^ The error bars indicate the 68 % confidence interval of 
the photon distribution fiP^. 
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